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ABSTRACT 

A  layer-ocean  model  was  developed  based  on  the  mixed 
layer  models  of  Kraus  and  Turner  (1957)  and  Denman  (1972) 
The  ocean  model  was  coupled  to  Pearson's  (1972)  time  vary- 
ing, symmetr ical- s tat ionary  hurricane  model  which  was 
based  on  a  model  proposed  by  Riehl  (1963) ,  and  in  which 
the  air-sea  interaction  was  specified  by  Cardone ' s  (1969) 
extension  of  Blackadar ' s  (1965)  two-layer,  baroclinic 
boundary- lay er  model.   Time-dependent  solutions  for  ocean 
mixed  layer  depth  and  temperature  were  obtained  in  res- 
ponse to  interaction  with  the  atmospheric  model.   Solu- 
tions indicated  that  the  interaction  between  entrainment 
mixing  and  upwelling  was  most  important  in  changing  mixed- 
layer  depth  and  temperature.   Radiational  effects,  inter- 
nal waves,  turbulent  scale  energy  dissipation,  and  large- 
scale  ocean  currents  were  not  included  in  the  present 
ocean  model . 
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I.    INTRODUCTION 


It  is  widely  accepted  that  the  source  region  for  the 
hurricane  heat  engine  is  the  underlying  ocean.   Although 
the  ocean  is  virtually  an  infinite  heat  source,  it  is 
generally  accepted  that  sea-surface  temperatures  in  excess 
of  26C  are  necessary  for  forming  and  maintaining  a  mature 
hurricane.   At  the  same  time  the  ocean  acts  as  a  sink  for 
the  hurricane's  momentum  which  is  imparted  to  the  ocean 
through  surface  stress.   It  is  the  balance  of  these,  two 
processes  which  was  the  basis  for  Riehl's  (1963)  steady- 
state  hurricane  model  in  which  the  tangential  wind  profile, 

outside  the  radius  of  maximum  winds,  was  specified  as 

1/2 

v_ r     =  constant . 
o 

Corgnati  (1971)  adopted  Riehl's  work  to  develop  a 
model  which  used  boundary-layer  predictive  equations  for 
moisture  and  temperature,  and  bulk  aerodynamic  equations 
for  the  fluxes  at  the  air-sea  interface.   Corgnati  in- 
cluded a  fairly  simple  ocean  model  which  described  the 
oceanic  response  to  convective  mixing  as  a  result  of  heat 
extraction;  however,  most  experiments  were  run  with  con- 
stant sea  surface  temperatures. 

Pearson  (1972)  modified  Corgnati's  model  by  using 
Cardone's  (1969)  extension  of  Blackadar's  (1965)  two- 
layer  near-neutral,  baroclinic  boundary -lay er  model  to 


calculate  the  fluxes  at  the  air-sea  interface.   Pearson's 
model  produced  steady-state  solutions  for  constant  sea- 
surface  temperatures,  and  exhibited  greater  sensitivity 
to  these  temperatures  than  did  the  Corgnati  model. 

It  was  felt  that  to  more  realistically  describe  the 
hurricane-ocean  interaction,  a  more  complex  ocean  model 
incorporating  wind  mixing,  convective  mixing,  and  advec- 
tion,  was  required.   It  was  the  object  of  this  study  to 
develop  such  a  model  and  couple  it  to  Pearson's  atmos- 
pheric model  through  the  boundary-layer  fluxes  of  heat 
and  momentum.   The  current  ocean  model  is  based  on  mixed- 
layer  models  of  Kraus  and  Turner  (1967)  ,  which  predicts 
the  mixed-layer  temperature  and  depth  in  response  to  wind 
mixing  and  convection,  and  Denman  (1972) ,  which  includes 
the  effect  of  upwelling.   It  was  hoped  that  the  ocean 
model  may  help  to  explain  which  of  the  mechanisms  is  of 
greatest  importance  in  producing  the  areas  of  low  sea- 
surface  temperatures  observed  in  the  wake  of  hurricanes. 
Investigators  have  attributed  the  low  temperatures  to 
various  mechanisms  with  little  agreement  as  to  which  is 
the  most  important.   For  example,  Jordan  (1954)  concluded 
that  mixing  was  the  more  important  factor.   Fischer  (1958) 
considered  upwelling  to  be  the  most  important.   Leipper 
(1967)  concluded  that  upwelling,  produced  by  surface 
divergence,  was  the  most  important  near  the  hurricane's 
center  and  that  mechanical  mixing  dominated  at  larger 
radii . 


The  present  ocean  model  was  run  using  each  of  the 
three  cooling  mechanisms  (convection,  wind  mixing,  and 
advection) ,  separately  and  in  various  combinations,  in 
an  attempt  to  study  their  relative  effects  upon  the 
ocean  upper  thermal  structure  and  subsequent  heat  content 
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II.    DESCRIPTION  OF  THE  MODEL 


A.   ATMOSPHERIC  VORTEX 

Riehl  (1963),  in  his  early  work,  assumed  conservation 
of  potential  vorticity  in  the  inflow  layer  and  a  bulk 
aerodynamic  equation  for  surface  stress  with  constant  drag 
coefficient  and  inflow  angle,  to  derive  the  relation 


1/2         ^   ^ 
v_r     =  constant 


(1) 


which  specifies  the  radial  profile  of  cyclonic  tangential 

winds  at  the  top  of  the  atmospheric  boundary  layer  as 

shown  in  Figure  (1) .   The  outer  limit  of  the  circulation 

described  bv  Riehl's  model  was  defined  by  the  radius  (r  ) 

o 

at  which  the  cyclonic  flow  becomes  anticyclonic  in  the 
outflow  layer.   Riehl  assumed  that  the  air  ascending  in 
the  eye-wall  region,  and  flowing  cyclonically  outward 
aloft,  conserves  absolute  angular  momentum.   With  this 
assumption,  one  obtains 


,2  ,1/2 

o      f    D  .   1 

1 


(2) 


where  f  is  the  Coriolis  parameter  and  the  i  subscripts 

denote  values  at  the  radius  of  maximum  wind  speed.   The 

present  atmospheric  model  assumed  that  the  radius  of 

maximum  wind  (r.)  was  coincident  with  the  eye-wall  radius 

1 

Since  the  eye  wall  is  a  preferred  region  for  deep  convec- 
tion and  latent  heat  release  in  observed  hurricanes, 
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those  boundary- lay er  air  parcels  with  the  greatest  heat 

and  moisture  content  in  the  model  were  assumed  to  be  at 

the  eye  wall.   Thus,  r.  was  taken  to  be  the  radius  at 

1 

which  the  equivalent  potential  temperature  (6  ) ,  averaged 
over  five  grid  increments  (Ar  =  3km) ,  was  a  maximum. 

Between  the  vortex  center  and  r.  the  winds  were  assumed 

1 

to  obey  v  r    =  constant,  which  agreed  with  data  from 

y 

Riehl  (1963),  as  well  as  Shea  and  Gray  (1972). 

From  observations,  Riehl  developed  the  relationship 


p'  =  -2.56  e 


(3) 


where  P'  is  the  pressure  departure  from  1005  mb  and  6  ' 

e 

is  the  equivalent  potential  temperature  departure  from 
350  K.   It  was  assumed  that  the  atmospheric  vortex  was 
in  a  state  of  quas i-equi libr ium  with  the  boundary  layer 
fluxes  of  heat. and  moisture,  which  determined  the  in- 
crease in  0  '  and  thereby  P'  of  the  inflowing  air. 
e 

Assuming  a  gradient  balance  between  the  wind  and 
pressure  field 

2 


1  9JP 

p  ar 


-^^  +  f  v„ 
r        9 


(4) 


and  using  the  relation 


V, 


2 

^e.   ^i 

1 


(5) 
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integration  of  (4)  from  r   to  r . ,  yields  the  following 
expression  for  the  maximum  wind  speed: 

0.5r 

v_   =  [  -  2  f  /r  .  r   -  r. 

6  .    r  -r .         1  o     1 
1      o   1 


+  /4f 2 (/r .r   -  r.)^  +  10^ 

1  o      1 


r  -r  .   (0 
o   1    e 

pr_ 


)  ] 


(6) 


For  the  symmetric  model,  the  time-dependent  predictive 
equations  for  potential  temperature  (6)  and  specific 
humidity  (q) ,  which  characterize  the  boundary-layer  pro- 
perties through  a  specified  depth  (AP),  are 


_9e 

9t 


9 (rv  9) 
r 

r  9r 


8  (0)9) 
8p 


<3Q. 


c  Ap 
p 


and 


(7) 


9q 

8t 


9 (rv  q) 

r 

r  9r 


9(ooq)  ^  ^^E 
9p      lAp 


(8) 


where  iQa)     ^^'^     ^^e^  ^^^    ^^^  oceanic  sensible  and  latent 

heat  fluxes,  respectively.   C   and  L  are  the  specific  heat 

and  latent  heat  of  vaporization.   Radial  (v  )  and  vertical 

r 

(CO)  components  of  motion  are  related  through  the  con- 
tinuity equation 


_9a)  _  _  r_ 

3p  ~     r  9r 


(9) 


The  surface  fluxes  Q  f  Q  /  T.  ,  and  the  integrated  radial 

S    E    9^ 

motion  (v  )  are  determined  from  the  boundary- lay er  model. 
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B.   BOUNDARY-LAYER  MODEL 

Following  Pearson  (1972),  the  boundary- lay er  model 
was  a  two-layer,  baroclinic  model  developed  by  Cardone 
(1969)  as  shown  in  Figure  (2) .   The  Cardone  model  of  the 
marine  boundary  layer  was  an  extension  of  the  near- 
neutral,  fixed- terrain  model  described  by  Blackadar  (1965)  . 

The  model  used  the  v_  profile  from  the  atmospheric  vortex 

U 

model,  together  with  a  computed  boundary- layer  thermal 
wind  and  made  use  of  Monin-Obuhkov  similarity  theory  and 
stationary  Ekman-layer  theory  to  calculate  the  friction 
velocity  u^,  the  stability  index  L^,  and  the  inflow  angle 
a.   These  are  then  used  to  calculate  a  vertical  mean  radial 


wind,  u,  within  the  surface  layer  from 

u  =  <  L   -  [Ln  ^  -  i){Lj] 
i  =  l  o 


/4 


(10) 


and  the  total  heat  flux  (H)  from 


H 


^s  "^  ^E  ^  [-u^C  pa^e]/KgL' 


(11) 


where 


1/2 

u^=  friction  velocity,   (T/p)     ,  where  T  is 

the  radial  stress  component, 

Z  =  6.84  X  10~^   +4.28  x  10~^  u^  -4.43  x  10~^, 
o   * 

^* 

surface  roughness  (m) , 

l|;(L^)=  wind  profile  function  [Cardone  (1969)], 

L  =  Z/L',  stability  index, 

Z.=  heights  in  the  surface  layer, 
1 
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K  =  von  Karman  constant, 

6  =  mean  potential  temperature  in  the 
boundary  layer,  and 

a, =  K, /K  ,  ratio  of  heat  transfer  coef f i- 
cient  to  turbulent  transfer 
coefficient  for  momentum. 

The  individual  Q   amd  Q   fluxes  were  specified  accord- 

s        E 

ing  to  the  Bowen  ratio  C  (0  -Q  ) /L  (q  -q  ) .   The  radial 

paw    V   a   w 

wind  within  the  spiral  layer  was  determined  through  in- 
tegration of  the  Ekman  solution  to  get 


/   V  dZ  =  /2    v^  sin  a 
h    ^  ^ 


___[eos(a+-  -Bh) 


/     37T     ,  .  , 

sm  (a+  —  -Bh)  ] 


(12) 


where 


B  =  (f/2u^Kh) 


1/2 


To  include  the  centripetal  acceleration,  the  Coriolis 

parameter  (f )  was  replaced  by  (f  +  ) .   The  values  of 

V   in  the  surface  layer  (Equation  10)  and  Ekman  layer 
(Equation  12)  were  then  used  in  determining  the  vertical 
motion  (Equation  9)  and  the  advection  of  heat  and 
moisture  [Equations  (7)  and  (8)]  in  the  prediction 
portion  of  the  model. 

C.    OCEAN  MODEL 

1 .   Description 

In  an  attempt  to  more  realistically  describe  the 
air-sea  interaction  associated  with  a  time-varying 
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hurricane  model,  it  was  decided  to  develop  an  ocean  model 
capable  of  producing  time-varying  solutions  of  mixed- 
layer  depth  and  temperature  in  response  to  atmospheric 
forcing.   Previously  the  hurricane  model  produced  time- 
varying  solutions  for  constant  sea-surface  temperatures 
until  a  steady  state  was  reached. 

The  present  ocean  model  was  initialized  with  sea- 
surface  temperature  (T  ) ,  mixed-layer  depth  (h  ) ,  below 

s  o 

layer  gradient  (9t/8z),  and  deep-layer  temperature  (T  ) 
as  shown  in  Figure  (3) .   By  allowing  the  ocean  thermal 
structure  to  change,  it  was  hoped  that  a  more  thorough 
understanding  of  the  air-sea  interaction  process  might 
be  possible.   The  present  model  is  designed  for  a  wind- 
dominated  region  where  there  is  negligible  downward  heat 
flux,  either  by  solar  radiation  or  sensible  and  latent 
heat  (condensation)  fluxes  from  the  atmosphere  to  the 
ocean.   In  the  absence  of  upwelling  this  precludes  any 
possibility  for  a  decreasing  mixed-layer  depth.   The 
atmospheric  forcing  of  the  ocean  model  was  by  surface 
kinetic  energy  input  (wind  stress)  and  upward  sensible 
and  latent  heat  (evaporation)  fluxes.   Both  wind  stress 
and  heat  flux  produced  changes  in  mixed-layer  depth  and 
temperature,  surface  stress  by  entrainment  mixing  and  the 
heat  fluxes  by  convective  mixing.   The  entrainment  mixing 
was  caused  by  surface  generated  turbulent  motions  at  the 
depth  of  the  mixed-layer  interface. 
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A  partitioning  procedure  was  used  to  budget  the 

total  atmospheric  stress  (t  )  between  surface  turbulence 

a 

production  and  generation  of  a  radial  mass  transport  ac- 
cording to  Ekman  theory.   The  effects  of  entrainment 
mixing  (turbulence)  were  assumed  to  be  redistributed 
uniformly  throughout  the  layer  by  turbulent  diffusion. 
A  radial  wind-driven  current  was  computed  from  the  mass 
transport  and,  by  employing  the  continuity  equation,  areas 
of  upwelling  and  downwelling  were  calculated.   Thus,  the 
ocean  model  incorporated  the  primary  mechanisms  believed 
to  produce  the  observed  ocean  cooling  in  the  wake  of 
hurricanes.   Internal  waves,  large-scale  currents,  and 
turbulent-scale  energy  dissipation  within  the  mixed  layer 
were  not  included.   At  six-hour  intervals  the  ocean  heat 
content  was  computed  and  compared  to  the  initial  values 
as  a  check  on  the  model's  solution  of  mixed-layer  depth 
and  temperature. 

The  time-dependent  mixed-layer  temperatures 
generated  by  the  ocean  model  produced  variations  in  the 
atmospheric  model's  v^  profile,  which  in  turn  produced 

V 

variations  in  the  atmospheric  forcing  (surface  stress  and 
heat  fluxes) .   This  tended  to  bring  the  coupled  models 
into  a  quasi-steady  state  condition  in  all  cases  except 
ones  in  which  advection  was  included. 

2 .   Mathematical  Development  of  Ocean  Model 

By  neglecting  molecular  heat  fluxes  and  viscous 
generation  of  heat,  Denman  (1972)  wrote  the  first  law 
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of  thermodynamics  (conservation  of  thermal  energy)  for 
an  incompressible  fluid  in  the  following  form: 


d_T 
dt 


P  c 
w  p 


(13) 


where 

—  Y  Z 
Q   =  Y^*^    /  specifies  the  heat  source  term,    (14) 

p   =  ocean  water  density, 
w 

y  =  average  extinction  coefficient, 
C   =  specific  heat  at  constant  pressure, 

P  JT  » 

R^  =  solar  radiation  incident  on  the  sea  surface, 
and 

Z   =  depth  in  the  ocean,  positive  downward. 

Substitution  of  (14)  into  (13)  gives  the  first  law  in 

the  form 


dT     „  -YZ 


(15) 


where 


R  =  R^/pC 


The  total  derivative  was  expanded  in  flux  form  and  (15) 
was  written  in  time-averaged  turbulent  form: 


8T    1  ^(^j-^T)     g  -  -     8(w'T')       -YZ 

■5T-  + ^ +  -^^ +  kt; =  YRe 

dt    r     dr         dZ  dZ 


(16) 


or 


dT 

—  =  YRe 


-YZ    8 (w'T' ) 


9z 


(17) 


where  u   =  vertical  average  radial  velocity  in  the  mixed 

layer,  and  T,  u   and  w  were  replaced  by  T  +  T',  u   +  u', 

r  r-         J  r     r 
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and  w  +  w',  respectively.   The  term  involving  (u'T')  in 
Equation  (16)  was  dropped  because  it  was  assumed  that 
there  was  no  systematic  correlation  between  u'  and  T' . 
Equation  (17)  is  the  heat  conservation  equation  where 


8(w'T')/8z  represents  the  local  divergence  of  the  turbu- 
lent vertical  heat  flux.   It  was  this  term  which  redis- 
tributed the  heat  exchanged  at  the  boundaries  uniformly 
throughout  the  homogeneous  mixed  layer.   At  the  ocean 


surface,  the  downward  turbulent  heat  flux,   (w'T')   ,  was 

o 

equal  to  the  nonradiative  heat  transfer  through  the 
ocean  surface,  i.e., 


(w'T' ) 


(Q„  +  Qo)/P  C 
^E    ^S    w  p 


(18) 


where  Q   and  Q   were  the  upward  latent  and  sensible  heat 
E        S 

fluxes,  respectively,  computed  in  the  boundary -layer 
model.   At  the  bottom  of  the  mixed  layer,  the  downward 
turbulent  heat  flux  was  equal  to  the  heat  transfer  through 
the  interface  due  to  entrainment  mixing,  i.e., 

dh 


where 


dh 


dt 


(w'T'),    =  (— ^  -  w)  (T^-  T,  ) 
h      dt  S    h 

o 


=  time  variation  of  mixed-layer  depth. 


(19) 


w  =  vertical  advection  velocity  at  depth  (h  ) , 

o 

T   =  mixed-layer  temperature,  and 

T   =  temperature  at  a  finite  depth  below  mixed- 
layer  depth. 


For  entrainment  mixing,  the  turbulent  flux  (w'T'),    was 

h 
o 

positive  which  corresponds  to  warm  water  (T'>  0)  being 
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forced  downward  (w'>  0)  and  cold  water  (T'<  0)  being 
lifted  upward  (w'<  0) .   This  process  was  accomplished 
against  buoyancy  forces  by  the  kinetic  energy  input  at 
the  surface  and  the  potential  to  kinetic  energy  conver- 
sion by  convective  overturning. 

Integrating  Equation  (17)  over  the  depth  of  the 

mixed  layer  (h  )  and  replacing  the  boundary  turbulent 
o 

heat  fluxes  by  Equations  (18)  and  (19),  yields  a  thermal 
energy  equation  for  the  mixed  layer, 


dT. 


-yh, 


dh 


^^h^  =  R(i-e   )-(Q^-^Qs^/P^%-(^  -)(^s-\^         ^2°^ 


By  expanding  the  total  derivatives  in  flux  form  and  ig- 
noring the  radiation  effects,   (R) ,   (20)  becomes  a  predic- 
tive equation  for  the  local  change  in  the  mixed-layer 

depth  (h  ) ,  i.e., 
o 


8h 
_( 

dt 


-h     9t^    ,  d (u  rT^)    3 (w  T^1 
o    S    _1     r   S  S 

(T^-  T,  )  ^dt       "*"  r    3r       "*"   3Z      ^ 


^E     S 


p  C  (T^-T,  ) 
w  p   S   h 


9h 
( 


+  w 


(21) 


To  develop  a  corresponding  equation  for  mixed-layer  tem- 
perature changes,  Denman  (1972)  integrated  Equation  (17) 
twice,  once  over  a  depth  z  and  once  over  the  depth  of  the 
mixed  layer  (h  ).   The  resulting  expression. 
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dT   h 
S  _o^ 

dt   2 


P  c 
w  p 


(w'T* ) 


dZ 


(22) 


represents  a  mechanical  energy  equation  for  the  mixed 

h    

layer.   The  term   /     (w'T')dZ  specifies  the  conversion 

o 
of  potential  energy  into  kinetic  energy  by  convection 

within  the  mixed  layer.   If  one  considers  the  mechanical 

energy  balance  expression 


W  +  G  -  D  =  0 


(23) 


where    G  =  kinetic  input  from  the  wind, 

D  =  dissipation  within  the  mixed  layer,  and  ■ 

h    

W  ="/    (w'T')dZ(potential  to  kinetic  energy 
o  conversion ) , 

substitution  of  Equation  (22)  into  (23)  yields 


dT^  h  ' 
S  _o_ 

dt    2 


P  c 
w  p 


h   -  (G-D) 
o 


(24) 


Expanding  the  total  derivative  of  temperature  in  flux 
form  and  solving  for  the  local  temperature  change  yields 


8t 


2 (G-D) 


h   p  c 
o   w  p 


1  ^^-r^^) 


dr 


d  (w  Tg) 
~8Z 


(25) 


which  is  a  mechanical  energy  equation  specifying  the 
change  in  mixed-layer  temperature  due  to  entrainment 
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mixing,  convection,  and  advection.   Equation  (25)  can 
be  used  to  eliminate  8t  /dt    from  the  heat  equation  (21) 


8h 
< 


2 (G-D) 


+ 


<2e-^  2s) 


h  (T  -T,  )    p  c  (T  -T,  ) 
o   s   h     w  p  s   h 


-  u 


9h 
< 

r  8r 


+  w 


(26) 


Equations  (25)  and  (26)  represent  a  system  of  two  equa- 
tions with  three  unknowns,  T  ,  h  ,  and  T,   (u   and  w  will 

so         h    r 

be  obtained  from  Ekman  theory  and  the  equation  of  con- 
tinuity) .   To  close  this  system  therefore,  a  third  equa- 
tion or  relation  was  required  which  specified  changes  in 
T,   as  being  dependent  upon  changes  in  h  .   The  parameter 

h  J      Jr  ir  3  Q  jr 

which  imposed  this  dependency  was  the  below-layer  gradient 
(8t/8z)  which  was  assumed  to  be  uniformly  affected  by  ver- 
tical velocities  and  thus  remained  constant.   Time  varia- 
tions of  T,  were  specified  by  the  expression 
h 


9t^  dh    .^ 

h      -    o,  ,9t 

■r— -  =  -  (w  -  TT— )  (tt^u 
dt  dt    dt  h 


(27) 


where 


w  =  vertical  advection  velocity 


The  solution  for  T,  ,  which  represented  the  lower  tempera- 
ture of  the  entrained  water,  proved  to  be  very  crucial 
in  the  mixed-layer  depth  and  temperature  solutions.   It 

was  the  temperature  difference  (T  -  T  )  which  imposed  the 

s    h 

primary  retarding  force  against  unlimited  deepening.   As 

time  progressed  the  difference  (T  -  T,  )  increased,  forming 

s    h 

an  actual  step  in  the  temperature  profile  at  the  base  of 
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the  mixed  layer.   Because  of  the  stable  lapse-rate 

(8t/8Z  <  0),  upwelling  (w<0)  contributed  to  a  lowering 

of  T   and,  hence,  increased  the  difference  (T  -  T,  )  . 
h  s    h 

Thus,  upwelling  had  the  effect  of  retarding  the  tendency 
for  mixed-layer  deepening.   At  the  same  time  upwelling 
(w<0)  enhanced  mixed- layer  temperature  changes  by  de- 
creasing the  layer  thickness  (Equation  26)  so  that  the 
energy  input  and  heat  fluxes  would  have  a  greater  effect 
(Equation  25).   Thus,  Equations  (25)   (26),  and  (27) 

formed  a  system  for  unknowns  h  ,  T  ,  and  T,  .   Space  and 

OS         h 

time  derivatives  in  (25)  and  (26)  were  approximated  by 
centered  differences,  except  for  the  initial  time  step 
which  used  a  forward  difference.   A  forward  step  was 
also  used  every  17th  time  step  (3  hours)  to  prevent  di- 
verging solutions.   Forward  time  differences  were  used 
in  (27).   The  difference  forms  of  (26)  and  (27)  were 

solved  together  for  h   and  T,  using  a  first  guess  depth 

oh 

and  temperature  (T  )  from  the  previous  time  step.   The 
first  guess  depth  in  (26)  and  (27)  was  increased  by 
0.05m  increments  until  the  difference  between  the  right 
and  left  hand  sides  of  the  difference  form  of  (26)  was 
less  than  or  equal  to  0.1m.   Using  the  new  mixed-layer 
depth,  mixed-layer  temperature  was  calculated  from  (25) . 

The  mechanical  energy  input  (G) ,  in  the  ocean 
predictive  equations  (25)  and  (26),  was  specified 
according  to  Kraus  and  Turner  (1965)  as 
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G  = 


gap. 


/a.  3/2 


(28) 


w 


where     g  =  gravitational  acceleration, 

1     w 

a  =  -t ,  coefficient  of  expansion, 

p   dT  tf  ' 

w 


air  density, 


p   =  sea  water  density, 
w 


u 


atmospheric  friction  velocity,  and 


T   =  atmospheric  surface  stress, 
o 

In  the  ocean  model,  surface  stress  (x  )  used  for  the 

o 

energy  input  term  (G)  depended  exponentially  upon  the 

ratio  of  the  mixed-layer  deoth  (h  )  and  the  roughness 

o 

length  (Z  ) : 
o 


T   =  T    EXP   (-^^7-) 
o     a        CZ 

o 


(29) 


where 


T   =  total  stress  available  from  boundary -layer 
model ,  and 

C   =  an  arbitrary  constant. 

h 
For  a  given  roughness  length  (Z  )  ,  the  function  exp  ("p"^ — ) 

o 
specified  the  percentage  of  total  stress  (T  ).  available 

for  turbulence  generation  as  a  function  of  mixed -layer 
depth  (Table  I) .   The  use  of  the  above  partitioning  pro- 
cedure was  an  attempt  to  budget  the  available  energy 
between  turbulence  generation  (G)  and  production  of  wind- 
driven  currents  in  the  ocean  model.   Experiments  were  run 
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with  various  values  of  C  and  compared  to  each  other  as 

well  as  to  a  run  in  which  all  the  surface  stress  was 

utilized  for  turbulence  generation.   The  dissipation,  D, 

was  set  to  zero  in  all  experiments. 

The  radial  velocity,  u  ,  which  is  needed  to  close 

r 

the  system,  was  calculated  in  accordance  with  Ekman  theory 
which  specifies  a  spiral  layer  to  a  depth  (Z  )  where  the 
current  direction  has  reversed  from  its  surface  direction 
and  the  speed  is  exp  (-7T)  of  the  surface  speed.   Below 
the  Ekman  layer  it  was  assumed  that  the  wind  had  no  direct 
influence.   According  to  Ekman  the  radial  component  of 
motion  within  the  spiral  layer  was  described  by 


u 


(t  -T  ) cos  a 
a   o 

P  (2f  iL_)V2 
w     p 
w 


EXP  (-BZ)  [cos  (-BZ) -  sin(-BZ)] 


(30) 


where    Z   =  Tr/B,  specifies  the  Ekman  depth. 


2y/p 


w 


f   =  coriolis  parameter, 

p   =  sea  water  density, 
w 

y   =  vertical  eddy  coefficient  of  viscosity, 

Z   =  depth,  and 

a   =  atmospheric  inflow  angle. 
Integration  of  (30)  over  a  specified  Ekman  depth  (Z  ) 
yields  an  expression  for  the  net  mass  transport  in  the 
radial  direction  for  the  Ekman  layer. 


M   =  I^/f 
r     6 


(31) 


25 


TABLE  I 


Stress  used  for  mixing 
mixed-layer  depth  as  a 
total  stress  (x  ) .   Z 
C  =10   and  5  x  ^10  .  ° 


(T  ) 
„  o 


vs 
of 


=  10  ^m 


Mixed  layer 
Depth  (m) 


as  % ,of  T 
°(C=10^     ^ 


T   as  %  of  T 
°(C  =  5x  10.^)^ 


5 
10 
20 
30 
40 
50 
60 
70 
80 
90 
100 


.95 
.90 
.82 
.  74 
.67 
.  61 
.  55 
.50 
.45 
.41 
.37 


.  90 
.82 
.67 

.  55 
.45 
.37 
.  30 
.  25 
.  20 
.  17 
.  14 
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where 


(t  -  T  ) cos  a,  tangential  atmospheric  surface 
a    o   ^ 

stress . 


This  represents  a  radial  mass  transport  per  unit  distance 

in  the  tangential  direction.   One  can  calculate  a  mean 

radial  advection  velocity  (u  )  for  the  mixed  layer  from 

r 

the  relation 


u  Z   =  M  /p 
re     r   w 


(32) 


The  mean  radial  (u  )  and  vertical  (w)  components  of  motion 

r 

are  related  through  the  continuity  equation 


8  (P   w) 
w 

9Z 


,  9  (p   u   r) 
—  w   r 

r    8r 


(33) 


Equation  (33)  was  integrated  over  depth  (Z  )  using  the 
boundary  condition  w(z=0)  =  0,  to  obtain  w  at  Z 


1  8     ,"  "^e, 

P   w  =  -  -  3—  (— z — ) 
w        r  dr    f 


(34) 


which  was  taken  to  be  w  at  h  .   The  probable  occurrence  of 

o 

a  return  flow  at  some  depth  below  the  mixed  layer  was  not 
included  in  the  present  ocean  model.   Thus,  the  only  effect 
on  the  below-layer  gradient  (9t/8z)  was  that  due  to  verti- 
cal velocities  and,  hence,  the  gradient  remained  constant. 
As  a  check  on  the  cooling  processes  (wind  mixing, 

convective  mixing,  and  advection)  ,  the  initial  ocean  heat 

2 
content  (cals/cm  )  for  a  depth  Z^  =  130m  was  calculated 

b 

based  on  the  initial  temperature  profile  (Figure  3) .   This 
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value  was  then  compared  to  the  heat  content  at  each  grid 
point,  computed  at  6-hour  intervals.   For  convective  mix- 
ing-only, with  the  kinetic  energy  input  (G)  and  the 
advection  velocities  (u   and  w)  equal  to  zero,  Equations 
(25)  and  (26)  predicted  excessive  mixed-layer  deepening 
and  cooling.   Also  the  final  grid  point  heat  contents  did 
not  reflect  just  the  heat  (H)  lost  to  the  atmosphere  due 

to  upward  heat  fluxes  ( Q   and  Q  ) .   Hence,  for  this  case 

E        S 

the  mixed-layer  depth  and  temperature  were  specified  by  a 
thermal  equilibrium  problem  using  the  following  set  of 
equations : 

T    -T     ' 

h    T       =    h     '      (T     '-    OAt/p    C    h')+    Ah       ^^    ^ 
OS  o  swpo  2 


h   =  h  •  +  Ah 
o     o 

T   =  T  '  +  4^  Ah 

s     s  8z 


(35) 


where 


h  '=  initial  mixed-layer  depth, 
o 

T  '=  initial  mixed-layer  temperature, 

8t/8z.  =  below-layer  gradient, 

Q    -    Q    "^    Q     >     total  upward  heat  flux, 
E     S 

p   =  sea  water  density, 
w 

At  =  time  step  (sec)  ,  and 

C   =  specific  heat 
P 

For  convective  mixing  only,  equations  (35)  were  solved  at 

each  time  step  for  the  change  in  mixed-layer  depth  (Ah) . 

With  this  calculated  depth  change,  the  new  layer  depth  (h  ) 

and  temperature  (T  )  were  specified. 

.   s 
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Figure  1.   Atmospheric  vortex  model 
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Figure  2.   Boundary  layer  model 
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Figure  3.   Ocean  model 
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Ill .    DISCUSSION  OF  RESULTS 

As  stated  previously,  the  object  of  this  study  was  to 
develop  an  ocean  model,  capable  of  producing  time-dependent 
solutions  for  mixed-layer  depth  and  temperature  in  response 
to  forcing  from  a  symmetrical,  stationary  hurricane  model. 
Previous  work  with  the  hurricane  model  was  done  with  con- 
stant sea-surface  temperatures  and,  consequently,  little 
was  learned  of  the  ocean  thermal  changes  or  of  the  corres- 
ponding changes  produced  in  the  hurricane  model. 

In  this  investigation  three  mechanisms  believed  to  con- 
tribute most  in  changing  the  ocean  thermal  structure  were 
isolated  and  studied.   These  were  entrainment  mixing  across 
a  stable  layer  (wind  stress),  convective  mixing  (latent 
and  sensible  heat  fluxes) ,  and  vertical  and  horizontal 
advection.   It  was  hoped  that  the  combined  time-dependent, 
atmosphere-ocean  model  would  help  determine  which  of  these 
mechanisms,  either  acting  alone  or  in  combination,  was  most 
important  in  producing  the  cold  sea-surface  temperatures 
observed  beneath  hurricanes.   To  accomplish  this,  a  series 
of  experiments  were  run  using  each  of  the  cooling  mechanisms 
alone  and  in  combination  with  each  other.   Figures  (4) 
through  (32)  depict  radial  profiles  of  various  ocean  para- 
meters for  the  different  experiments.   For  these  runs  the 
ocean  model  was  initialized  with  mixed-layer  depth  of 
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30  meters,  mixed-layer  temperature  of  30C,  below-layer 

gradient  of  lOC/lOOm,  and  deep-layer  temperature  of  20C. 

The  lower  boundary  of  the  model  was  fixed  at  130  meters 

(Figure  3) ,  below  which  there  were  assumed  to  be  no 

thermal  effects.   Other  initial  temperature  profiles  can 

be  used;  however,  the  above  values  were  believed  to  be 

fairly  representative  of  the  tropical  ocean  regions.   The 

atmospheric  model  was  initialized  with  air  temperature  of 

29C,  radius  of  maximum  wind  (r.  =  27km) ,  and  maximum  wind 

1 

(v^   =  30m/sec) .   In  runs  with  various  constant  sea-surface 
u  . 

1 

temperatures,  Pearson  (1972)  found  a  steady-state  v..   of 

U  . 

1 

60m/sec  for  30C. 

The  value  of  the  constant  (C)  in  the  wind  stress  parti- 

4 
tioning  equation  (29)  was  arbitrary  fixed  at  10   in  all 

the  initial  runs  except  the  advect ion-only  case.   For  this 

run  the  partitioning  process  was  not  used;  consequently, 

all  the  available  stress  (T  )  was  used  to  produce  radial 

a 

wind-driven  currents  (u  )  and  vertical  velocities  (w)  . 

r 

Since  the  partitioning  process  was  an  important  assumption 
in  the  model,  an  entr ainmen t-convec t ion-advec t ion  experi- 
ment was  run  with  a  value  of  5  x  10   for  the  constant  (C) 
in  Equation  (29) .   With  a  smaller  value  of  C,  Equation  (29) 
partitioned  more  atmospheric  stress  to  producing  ocean 
currents  and  less  to  entrainment  mixing  at  all  mixed-layer 
depths  (Table  I).   An  en trainment-convec t ion  run  was  also 
made  in  which  the  partitioning  process  was  not  used,  hence 
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all  available  stress  from  the  boundary- lay er  model  was 
used  in  the  kinetic  energy  term  (G)  in  the  mixed-layer 
depth  and  temperature  predictive  Equations  [(25)  and  (26)] 

A.  CONVECTIVE  MIXING  ONLY 

For  the  convective  mixing-only  case,  Equations  (35) 
were  solved  for  the  mixed-layer  depth  (Figure  4)  and  tem- 
perature (Figure  5) .   The  atmospheric  forcing  was  by 
latent  and  sensible  heat  fluxes  from  the  ocean  surface  to 
the  atmosphere  as  computed  in  the  boundary -lay er  model. 
Maximum  cooling  (0.5C)  and  deepening  (5  meters)  of  the 
mixed  layer  at  18  hours  occurred  under  the  region  of  maxi- 
mum wind  (r.)  as  expected  from  Equation  (11)  which  relates 
total  oceanic  heat  loss  (H)  to  the  cube  of  the  friction 
velocity  (u^ ) .   Oceanic  heat  loss  computed  from  initial 
and  final  temperature-depth  profiles  (Figure  6)  in  the 
ocean  model  compared  quite  favorably  with  the  total  heat 
extraction  (H)  computed  in  the  boundary- lay er  model 
(Figure  7) .   The  amount  of  ocean  cooling  due  to  convection 
was  comparable  to  Jordan's  (1964)  estimates. 

B.  ENTRAINMENT  MIXING 

Radial  profiles  of  mixed-layer  depth  and  temperature 
due  only  to  entrainment  mixing  (Figures  8  and  9)  were  ob- 
tained by  iterative  solutions  of  Equations  (25)  ,   (26)  , 

and  (27).   There  was  no  ocean  current  generation  in  this 

4 
run,  and  the  partitioned  stress  (29) ,  with  C  =  10  ,  was 
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Figure  4.   Radial  profile  of  mixed-layer  depth  at 
18  hours  due  to  convective  mixing. 
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temperature . 
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Figure  6.   Predicted  18-hour  vertical  temperature 
profile  at  30  km  for  convective  mixing. 
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ocean  heat  loss  and  accumulated  (Q  +0  ) 

E   S 
from  boundary  layer  model. 
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Figure  8.   Radia]  profiles  of  mixed-layer  depth  at 
6  and  18  hours  for  entrainment  mixing. 
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the  forcing  which  produced  entrainment  mixing  through  sur- 
face turbulent  generation  (G).   Maximum  cooling  (0.5C) 
and  deepening  (23  meters)  at  six  hours  and  1 . 5C  and  45 
meters  at  18  hours  again  occurred  beneath  the  region  of 
maximum  kinetic  energy  input  (28)  as  expected  from  (25) 
and  (25)  .   The  increased  cooling  and  deepening  v;ith  time 
at  larger  radial  distances  was  due  to  the  outv/ard  migra- 
tion of  the  hurricane  eye  wall  (r.)  in  response  to  colder 

1 

sea-surface  temperatures  near  the  vortex  center. 

Since  entrainment  mixing  produced  only  vertical  heat 
redistribution,  the  heat  content  of  the  ocean  beneath  the 
atmospheric  vortex  should  be  conserved  except  for  numeri- 
cal errors.   Figure  (10)  depicts  the  radial  distribution 
of  oceanic  heat  loss  at  18  hours  for  the  en tr a inment-only 

experiment.   The  area  averaged  heat  loss  for  the  solution 

2 
domain  (Table  II)  was  137  cals/cm  .   Considering  the 

sensitivity  of  heat  content  to  changes  in  temperature  and 

depth  and  since  the  depth  and  temperature  solutions  were 

computed  numerically,  and  are  subject  to  truncation  errors, 

the  values  in  Figure  (10)  are  considered  within  acceptable 

limits.   For  example,  in  the  region  of  maximum  mixed-layer 

depth  (75m),  a  temperature  difference  of  O.IC  will  give 

2 

a  heat  contant  change  of  750  cals/cm  .   The  maximum  error 

2 
in  Figure  (10)  is  800  cals/cm  . 

Figure  (11)  depicts  ocean  temperature  profiles  at  6, 

12,  and  18  hours  as  computed  in  the  model.   Note  the  strong 
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Predicted  18-hour  radial  profile  of 
ocean  heat  loss  for  entrainment  mixing 
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Figure  11 


Predicted  6-,  12-,  and  18-hour  vertical 
temperature  profiles  at  30  km  for 
entrainment  mixing. 
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gradient  which  develops  below  the  mixed  layer.   As  men- 
tioned in  a  previous  section,  the  temperature  difference 

(T  -  T,  )  is  the  primary  retarding  force  which  balances 
s    h 

the  atmospheric  forcing.   The  below-layer  gradient  {t: — ) 

o  Z 

is  an  important  factor  in  determining  the  depth  at  which 
this  balance  will  occur.   For  a  given  energy  input  (G) , 
one  expects  deeper  mixing  for  a  weak  gradient  and 

shallower  mixing  for  a  strong  gradient. 

/ 

C.    ADVECTION-ONLY  CASE 

For  the  advect ion-only  case  (G=0,  and  Q    -^    Q    =0)  ,  mixed" 

E    S 

layer  depth  and  temperature  changes  were  determined  by  the 

advection  terms  in  the  predictive  equations  [(25)  and  (26)] 

The  partitioning  process  (29)  was  not  used,  hence,  all 

available  atmospheric  stress  (T  )  was  used  to  produce  a 

a 

radial  mass  transport  (32)  and  a  mean  radial  wind-driven 

current  (33) .   The  rapid  increase  in  u   (Figure  12)  from 

r 

the  vortex  center  to  the  radius  of  maximum  wind  (r.)  pro- 

1 

duced  a  region  of  horizontal  divergence  and,  hence,  enor- 
mously strong  upwelling  (Figure  13) .   Just  beyond  radius 
r.  a  small  region  of  convergence  and,  hence,  downwelling 

occurred.   However,  the  gradual  decrease  in  u   was 

r 

balanced  by  the  increased  radius  to  maintain  a  fairly  con- 
stant total  mass  transport.   Consequently,  areas  of  signi- 
ficant downwelling  were  not  evident  within  the  solution 
domain  (300  km) . 

The  effects  of  strong  upwelling  were  quite  evident  in 
the  mixed-layer  depth  and  temperature  profiles  (Figures 
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14  and  15) .   Within  six  hours  the  mixed  layer,  near  the 
vortex  center,  was  eliminated  and  the  sea-surface  tempera- 
ture decreased  until  the  cold  deep  ocean  water  (20C)  was 
upwelled  to  the  surface.   At  18  hours  the  radius  of  maxi- 
mum wind  (r. )  had  moved  outward  in  response  to  cold  sea- 
1 

surface  temperatures  at  the  vortex  center  (Figure  16) . 
The  region  of  maximum  upwelling  likewise  moved  outward, 
leaving  a  broad  region  of  20C  sea-surface  temperatures 
between  15  and  45  km  and  zero  mixed-layer  depth  to  a 
radius  of  95  km.   The  dramatic  effects  of  advection  were 
also  evidenced  by  significant  ocean  heat  content  changes 
(Figure  17) .   The  positive  values  correspond  to  heat  loss 
regions  and  the  negative  values  to  heat  gain  regions. 
Vertical  and  radial  velocities  produced  a  redistribution 
of  heat  by  changing  the  ocean  temperature  profile  (Figure 
18) .   Because  of  the  ocean's  response  to  the  atmospheric 
stress  profile,  radial  currents  (u  )  transported  heat 
across  the  model's  open  boundary  at  300  km.   Consequently, 
there's  a  significant  averaged  heat  loss  (Table  II)  in 
the  ocean  model  for  this  experiment. 

D.    ENTRAINMENT  AND  CONVECTIVE  MIXING 

This  experiment  was  similar  to  the  entr a inment-on ly 

run  except  heat  flux  iQ    +    Q    )     was  included  in  (25)  and 

E    S 

(26).   For  the  initial  run,  stress  partitioning  (29)  was 

4 
utilized  with  C  =  10  .   Radial  profiles  of  mixed-layer 

depth  (Figure  19)  and  temperature  (Figure  20)  closely 
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Figure  16.   Same  as  Figure  12  except  tangential 
component  of  the  wind. 
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Figure  17.   Same  as  Figure  12  except  ocean  heat  loss. 
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Figure  18  . 


Predicted  6-  and  12-hour  vertical 
temperature  profiles  at  30  km  for 
advection  only. 
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resemble  the  corresponding  profiles  for  en tr a inmen t-only 
(Figures  8  and  9).   The  inclusion  of  the  heat  flux  term 
produced  a  slight  increase  in  the  amount  of  deepening  and 
cooling  compared  to  entr ainment-on ly ;  however,  the  amount 
was  less  than  2  meters  and  0.5C  at  18  hours.   One  may 
conclude  that  entrainment  mixing  dominates  when  compared 
to  convective  mixing  for  the  initial  conditions  and  atmos- 
pheric forcing  utilized  in  the  model. 

For  this  run  the  oceanic  heat  content  should  not  be 
conserved  because  of  the  latent  and  sensible  heat  exchange 
with  the  atmosphere.   Figure  (21)  depicts  radial  profiles 
of  oceanic  heat  loss,  calculated  in  the  ocean  model,  and 
heat  extraction  (H)  from  the  boundary-layer  model.   The 
difference  between  the  profiles  represents  the  loss  attri- 
buted to  inaccuracies  in  the  numerical  solution  of  (25) , 
(26) ,  and  (27) ,  and  unfortunately  is  as  large  as  H  itself. 
The  vertical  ocean  temperature  profiles  (Figure  22)  were 
similar  to  those  for  entr a inment-on ly  (Figure  11).   For 
comparison.  Figures  (23)  and  (24)  show  radial  profiles  of 
mixed-layer  depth  and  temperature,  respectively,  for  an 
entrainment-convec t ion  run  in  which  stress  partitioning 

was  not  used.   All  available  stress  (t  )  from  the  boundary' 

a 

layer  model  was  used  in  the  kinetic  energy  term  (G)  of 
(25)  and  (26)  to  produce  changes  in  mixed-layer  depth  and 
temperature.   The  increased  atmospheric  forcing  was 
evidenced  by  an  increased  maximum  mixed-layer  depth  (12  m 
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Figure  21.   Predicted  18-hour  radial  profiles  of 

ocean  heat  loss  and  accumulated  (Q  +Q„) 
from  boundary- layer  model. 
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Figure  22.   Same  as  Figure  11  except  for  entrainment 
and  convective  mixing- 
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Figure  23.   Radial  profile   of  mixed-layer  depth 
at  18  hours  for  entrainment  and 
convective  mixing  (no  partitioning). 
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deeper)  and  a  lower  minimum  mixed-layer  temperature  { . 5C 
colder)  than  for  the  partitioned  run.   In  both  of  these 
ent ra inment- convec t ion  runs,  as  well  as  the  entrainmcnt- 
only  run,  the  amount  of  sea-surface  cooling  near  the  vor- 
tex center  was  sufficient  to  cause  the  radius  of  maximum 

wind  (r.)  to  migrate  outward.   The  intensification  of  the 
1 

atmospheric  forcing  at  larger  radii  produced  a  correspond- 
ing increased  deepening  and  cooling  of  the  mixed  layer  at 
larger  radii.   Since  there  was  no  upwelling  in  these  runs, 
the  radial  profiles  of  mixed-layer  depth  and  temperature 
were  not  modified,  and  the  radius  of  maximum  cooling  and 
deepening  remained  at  20  km. 

In  the  next  set  of  experiments  advection  will  be  in- 
cluded and  the  effect  of  variations  in  the  partitioning 
constant  (C)  will  be  examined. 

E.    ENTRAINMENT,  CONVECTION,  AND  ADVECTION 

This  run  incorporated  all  the  mechanisms  believed  to 
be  of  importance  in  cooling  the  ocean  beneath  a  hurricane. 
The  total  available  atmospheric  stress  (x  )  was  partitioned 

Si 

(29)  into  that  which  produced  surface  turbulence  (T  )  and 

that  which  generated  wind-driven  currents  (T  -T  ) .   For 

a   o 

the  initial  run  the  partitioning  constant  (C)  in  Equation 

4 
(29)  was  set  equal  to  10  .   This  gave  a  stress  vs.  depth 

profile  as  depicted  in  Table  I.   In  addition  to  the  kinetic 

energy  and  flux  terms,  advective  terms  were  included  in  the 

solution  for  mixed-layer  depth  (26)  and  temperature  (25) . 
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The  radial  advection  velocity  (u  )  was  specified  according 

r 

to  (33)  and  areas  of  vertical  velocities  (w)  were  calculated 

from  (34).   Radial  profiles  of  u   (Figure  25)  and  w  (Figure 

25)  were  similar  in  nature  to  those  for  advection  only 

(Figures  12  and  13) .   However,  since  stress  partitioning  was 

not  used  for  the  advect ion-only  experiment,  the  magnitudes 

of  u   and  w  were  smaller  in  the  current  experiment. 
r 

As  in  previous  runs,  the  hurricane  eye  wall  moved  out- 
ward in  response  to  cold  sea-surface  temperatures  near  the 
vortex  center  (Figure  15) .   The  upwelling  region  (Figure  25) 
likewise  moved  outward  in  response  to  the  movement  of  maxi- 
mum surface  stress.   The  changes  in  mixed-layer  depth  and 
temperature  produced  by  entrainment  and  convective  mixing 
(Figures  19  and  20)  were  significantly  modified  by  the  out- 
ward migration  of  the  upwelled  region.   For  the  entrainment- 
convection  mixing  case  the  region  of  maximum  mixed-layer 
deepening  was  coincident  with  the  region  of  maximum  cooling. 
However,  with  the  addition  of  advection,  the  maximum  deepen- 
ing (Figure  27)  occurred  at  radii  just  outside  the  radius 
of  maximum  wind  (r.),  where  upwelling  was  weak  and  the  atmos- 
pheric forcing  (G  and  Q  +  Q  )  was  fairly  large.   At  smaller 

F     S 

radii  the  upwelling  dominated  the  effect  of  mixing,  leaving 
very  shallow  mixed-layer  depths.   The  region  of  maximum 
ocean  cooling  (Figure  28),  which  remained  in  the  vicinity 
of  15  km,  represented  an  interaction  between  mixing  and  ad- 
vection.  By  decreasing  the  layer  thickness,  upwelling 
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Figure  27.   Same  as  Figure  25  except  mixed-layer 
depth . 
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enhanced  the  effects  of  atmospheric  forcing  on  the  mixed- 
layer  cooling  (25).   Stress  partitioning  increased  the 
percentage  of  total  stress  used  in  the  kinetic  energy  term 
(G)  as  the  mixed-layer  depth  decreased.   There  was  a  corres- 
ponding decrease  in  the  stress  used  for  currents  (upwelling) 
which  tended  to  keep  the  mixed-layer  depth  from  becoming 
zero.   Since  ultimately  it  was  the  value  of  the  constant 
(C)  which  determined  whether  mixing  or  upwelling  would  have 
the  greatest  influence  on  the  mixed-layer  depth,  the  entire 
formulation  of  stress  partitioning  deserves  considerable 
future  study. 

Significant  ocean  heat  loss  (Figure  29)  occurred  inside 
30  km  where  upwelling  had  the  greatest  effect  on  the  verti- 
cal temperature  profile  (Figure  30) .   At  larger  radii,  where 
entrainment  and  convection  dominated,  ocean  heat  content 
changes  reflected  the  heat  extraction  (H)  by  the  boundary- 
layer  model.   The  area  averaged  heat  loss  (Table  II)  for 

2 
the  solution  domain  was  2300  cals/cm   at  18  hours.   As  in 

the  advection-only  case,  this  loss  reflects  the  transport 

of  heat  across  the  model's  open  boundary  at  300  km. 

To  test  the  influence  of  stress  partitioning  (29)  ,  a 

3 

second  experiment  was  run  using  a  value  of  C  =  5  x  10  . 

With  this  smaller  value,  the  percentage  of  total  stress 

(t  )  used  for  entrainment  mixing  was  smaller  and  decreased 
a 

more  rapidly  with  mixed-layer  depth  and,  consequently,  the 

percentage  of  stress  used  for  current  generation  (T  -T  ) 

a   o 
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Figure  29.   Predicted  6-  and  18-hour  radial  profiles 
of  ocean  heat  loss  for  entrainment, 
convection,  and  advection. 
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Figure  30.   Same  as  Figure  11  except  entrainment, 
convection,  and  advection. 


59 


was  larger  and  increased  more  rapidly  with  mixed-layer 

depth  (Table  I).   Comparing  radial  profiles  of  mixed- 

4 
layer  depth  and  temperature  for  C    -    10        [Figures  (27)  and 

3 
(28)]  with  those  for  C  =  5  x  10   [Figures  (31)  and  (32)], 

it  is  apparent  that  advection  by  horizontal  and  vertical 

velocities  had  a  greater  influence  in  the  latter  case. 

[Note  similarity  between  Figures  (31)  and  (32)  and  those 

3 

for  advection  only  Figures  (14)  and  (15)] .   For  C  =  5  x  10 

upwelling  has  eliminated  the  mixed  layer  in  a  small  region 
near  the  vortex  center  and  the  combined  effects  of  upwell- 
ing and  mixing  have  cooled  the  sea-surface  temperature  to 

20C.   At  larger  radii  where  mixing  dominates,  there  was 

4 
greater  deepening  and  cooling  with  C  =  10   than  with 

3 

C  =  5  x  10  . 

Table  II  summarizes  some  significant  atmospheric  and 

ocean  variables.   The  air-sea  interaction  was  evident  in 

all  cases  where  there  was  pronounced  cooling  of  the  ocean 

surface  near  the  vortex  center.   The  lower  sea-surface 

temperatures  reduced  the  air-sea  temperature  difference, 

hence  the  heat  fluxes  from  the  ocean  to  the  boundary  layer 

were  reduced.   The  tendency  for  the  model's  eye  wall  to 

seek  a  region  of  maximum  equivalent  potential  temperature 

(6  )  caused  the  radius  of  maximum  wind  (r.)  to  move  outward 
e  1 

as  the  inner  surface  was  cooled  (Figure  16) .   This  outward 

migration  maintained  a  sufficient  equivalent  potential 

temperature  difference  (G   -  9   )  and,  consequently,  the 

e  .     e 
1     o 
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Figure  31.   Same  as  Figure  27  except  C=5xl0' 
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Figure  32.   Same  as  Figure  28  except  C=5xlO' 
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maximum  wind  (v„  )  remained  fairly  constant.   Because  of 

0  . 

1/2        ■"" 
the  V  r     =  constant  profile,  the  outv/ard  migration  of 
u 

r.  produced  an  unrealistic  intensification  of  the  circula- 
1 

tion  at  large  radii.   This  in  turn  maintained  a  constant 
total  ocean  mass  transport  with  radius  which  resulted  in 
significant  amounts  of  heat  being  transported  out  the 
model's  300  km  boundary  to  a  region  where  downwelling  would 

presumably  have  occurred.   In  order  to  properly  model  the 

1/2 

movement  of  the  eye  wall  the  v«r     =  constant  profile  must 

D 

be  modified  to  prevent  this  unrealistic  intensification  at 
large  radii . 
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TABLE  II 

Significant  atmospheric  and  oceanic 
variables  for  different  experiments 
at  18  hours . 


max  min 

h  T 
o     s 
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Area  Aver- 
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»    2 
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4 
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IV.    SUMMARY  AND  CONCLUSIONS 

The  results  of  this  study  indicate  the  feasibility  of 
coupling  an  ocean  model  to  an  atmospheric  vortex  model  to 
simulate  time-dependent  air-sea  interaction.   Since  the 
model  was  stationary,  and  the  effects  of  large-scale  ocean 
current,  turbulent- seal e  energy  dissipation,  radiation, 
and  internal  waves  were  not  included,  the  magnitude  of  the 
ocean  thermal  changes  may  not  be  representative.   However, 
the  relative  changes  for  the  various  cooling  mechanisms 
(convective  mixing,  entrainment  mixing,  and  advection)  were 
apparent. 

First,  convective  mixing  had  insignificant  effects  on 
an  ocean  thermal  structure  initialized  with  30  meter  mixed- 
layer  depth,  30C  mixed-layer  temperature,  and  below-layer 
gradient  of  lOC/lOOm.   Convective  cooling  would  have  been 
more  pronounced  for  shallower  layer  depths  and  stronger 
below-layer  gradients;  however,  these  were  not  tested  in 
the  model. 

Second,  entrainment  (wind)  mixing  produced  fairly  sub- 
stantial mixed-layer  deepening  and  cooling  in  the  region 
of  maximum  wind  stress.   The  radius  of  maximum  wind  res- 
ponded to  the  colder  sea-surface  temperatures  by  migrating 
outward.   The  effect  of  entrainment  mixing,  likewise, 
moved  outward  producing  deeper  layer  depths  and  colder  sea- 
surface  temperatures  at  larger  radii.   The  depth  and 
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temperature  changes  due  to  entrainment  and  convection  were 
comparable  to  those  for  entrainment  only  indicating  en- 
trainment to  be  the  dominate  process. 

Third,  the  effects  of  the  wind-driven  circulation, 
when  acting  alone,  produced  the  most  pronounced  changes 
in  mixed-layer  depth  and  temperature.  ■  Upwelling  inside 
the  radius  of  maximum  wind  displaced  the  relatively  warm 
mixed-layer  water  by  cold  deep-ocean  water  advected  from 
below.   The  warm  sur face- lay er  water  was  advected  outward 
by  radial  currents.   Since  no  significant  downwelling 
occurred  within  the  solution  domain,  heat  was  lost  through 
the  model's  outer  boundary. 

In  reality  the  various  mechanisms  studied  do  not  act 
separately  but  rather  interact  in  some  complex  fashion  to 
produce  the  thermal  changes  which  observations  have  shown. 
It  appeared  to  be  a  question  between  entrainment  mixing 
and  advection  as  to  which  has  the  dominate  role  in  cooling 
sea-surface  temperatures  beneath  a  hurricane.   Advection 
definitely  produced  significant  cooling  inside  the  radius 
of  maximum  wind  but  to  think  of  a  wind-driven  circulation 
with  no  surface  turbulence  seems  unrealistic.   It  appeared 
that  the  interaction  of  the  two  processes,  with  a  budget- 
ing of  the  total  atmospheric  stress  between  turbulence 
generation  and  current  generation,  was  the  dominate 
mechanism  for  sea-surface  cooling. 
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